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A magnetic flux line lattice (FLL) or vortex lattice (VL) 
with lattice constant aooc ^ ieU was predicted jlj] in the 
entire mixed state of a type-II superconductor within a 
mean field model. The advent of high temperature su- 
perconductors (HTSC) focused attention on the influence 
of thermal fluctuations on VL and led to the prediction 
H of a Vortex Liquid State with an unusual reentrant 
phase boundary in the magnetic field-temperature (H, T) 
plane such that at a fixed T, the melting phase boundary 
(H m (T)) is encountered twice upon increasing H ||. A 
dilute vortex liquid phase is expected to form above the 
lower critical field H c i, where ao > A, the magnetic pene- 
tration depth A gives a measure of the range of repulsive 
interaction that governs the FLL. A highly dense vortex 
liquid phase is expected below the upper critical field, 
where ao > 2£; £ is the coherence length, i.e., the radius 
of the vortex core. The two branches comprising reen- 
trant melting phase boundary join around the so-called 
"nose" temperature above which the vortex solid phase 
is thermally melted at all field values. The upper branch 
of the melting curve is given approximately by eqn. 5.5 
ofRef.3, 



(1) 



where Cl is the Lindemann parameter, Gi is the Ginzburg 
number (= (l/2)(k B T c / H^e/STr) 1 / 2 ) and (3 m w 5.6 
. Gi is much larger in the Cuprate HTSC systems as 
compared to that in low T c alloys, which facilitates the 
observation of the dense vortex liquid phase in HTSC 



systems jj) . However an experimental observation of the 
dilute vortex liquid phase, has remained elusive. On the 
lower branch, the melting phase boundary is governed by 
eqn. 5.19 of Ref.3, 

-e ( ^=/U%^)(§) 2 (l-£) (2) 

where a o =[2/3 1 / 2 (0 o /B)] 1/2 and /3 m '=0.5. Solving 
eqn. (2) gives B m (T^O) = H cl (0) pn(T/T c )] " 2 . The large 
value of k in the Cuprates presumably restricts the di- 
lute vortex liquid phase to very small induction values, 
making its experimental observation a challenging task. 

A further difficulty in the observation of the vortex liq- 
uid phases arises due to the ubiquitous quenched disorder 
(i.e., pinning centers) in all real systems which adversely 
affects the stability of the ordered vortex solid phase. The 
effect of disorder is expected to be especially strong for 
dilute vortex arrays and may mask the observation of an 
underlying dilute liquid - dilute solid melting transition 

In the context of quenched disorder, the phenomenon 
of the peak effect (PE), an anomalous increase in the 
critical current, is continuing to receive a great deal of 
attention as it marks the competition between interac- 
tion and disorder (i.e., elastic energy versus pinning en- 
ergy) ||-[l|. In a collective pinning description |Q of 
superconductors, J c relates inversely to the volume V c of 
Larkin domain ( J r cx i- ) over which FLL is correlated. 

The anomalous peaking in J c indicates a rapid decrease 
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in V c due to softening of elastic moduli of FLL at the 
incipient melting transition |6| JTll , [l2] , [l^] , thus making the 
PE yet another signature of the loss of order of the vortex 
array. The identification of PE as a "phase boundary" 
remains among the more complex and controversial ques- 
tions in type II superconductors in general frll^]. 

Ghosh et al |10| showed that in a weakly pinned vor- 
tex lattice of hexagonal 2H-NbSe 2 (T c w7.1 K), the PE 
curve, which is the locus of peak temperatures T p (H) 
(determined from the magnetic shielding response, i.e., 
ac susceptibility \' measurements), bore a striking re- 
semblance to the theoretically proposed reentrant [^]]|] 
melting phase boundary T m (H). The turnaround of reen- 
trant T p (H) curve has been reported to occur when 
a (of 4000 A at H « 150 Oe) «A C (in 2H-NbSe 2 ) , and 
it is followed by a rapid broadening of PE. The PE ul- 
timately became undetectable below 30 Oe jio). Since 
PE is directly related to pinning, its variation with vary- 
ing pinning strength is crucial in providing an under- 
standing of the complex interplay between thermal and 
quenched disorder in destabilizing the ordered phase. In 
this letter, we report new results on the reentrant nature 
of PE curve and the effect of quenched random disor- 
der on it via electrical transport measurements and an 
equilibrium dc magnetization experiment on a specimen 
of 2H-NbSe2 utilized by Ghosh et al |Io| and magnetic 
shielding (%') response studies on two different crystals 
of the same compound with qualitatively different lev- 
els of quenched random disorder. The conventional dc 
electrical transport results not only fortify the earlier 
claim |l0| of reentrant nature of the T p (H) curve but 
also corroborate the previous finding that PE broad- 
ens across fields where turnaround in T p (H) takes place. 
The magnetization hysteresis data also confirm the dou- 
ble crossover of T p (H) curve in an isothermal scan near 
the turnaround region. We find that though the upper 
(higher H) branch of T p (H) curve is somewhat robust, the 
lower reentrant branch is strongly influenced by disorder. 
Decreasing pinning strength results in the 'nose' region 
being pushed down to fields lower than that reported 
earlier whereas increasing pinning could make the 
'nose' feature completely disappear. These findings pro- 
vide quantitative estimate on how increasing disorder 
shrinks the (H,T) space over which vortex array retains 
spatial order. 

2H-NbSe2 has an appreciable value («3 x 10 -4 ) of Gi, 
which lies between those of HTSC and conventional su- 
perconducting alloys [0] . The weak pinning character- 
istic of single crystals of this compound [J c /Jo^l0~ 5 to 
10 -7 , where J c is the critical current density and Jo (~ 
10 8 A/cm 2 ) is the theoretical depairing current density] 
facilitates the formation of VL in them |LlJ and easy vari- 
ability of spatial correlation lengths in its specimen make 
them convenient test beds to observe phenomena which 
are a consequence of competition and interplay between 
interaction, thermal fluctuations and quenched random 
disorder. Three single crystals with increasing quenched 
disorder used in the present study comprise platelets be- 



longing to the same batches as the specimen X used by 
Higgins and Bhattacharya [O], as the specimen Y used 
by Ghosh et al |10| and as the specimen Z used by Hen- 
derson et al |13|. The dc resistance data to ascertain 
PE temperatures in specimen Y were recorded following 
Ref. 11, the isothermal magnetization hysteresis data to 
identify the manifestation of PE in a similar sample were 
obtained on a Quantum Design Squid Inc. magnetome- 
ter following the half scan technique of Ravikumar et al 
]l5t and ac susceptibility data in all three samples were 
recorded using a high sensitivity ac susceptometer p6| . 
The essential new findings are summarized in Figs. 1 
to 3. All data presented are for H||c, however, PE is ob- 
served in all orientations in 2H-NbSe2 and the anisotropic 
Ginzburg-Landau description applies to this system [|ll| . 

Fig.l shows the temperature variation of resistance, R 
vs reduced temperature (t=T/T c (0)), in the crystal Y at 
some of the low fields (120 Oe to 380 Oc) to locate the 
signature of the reentrant behavior in T p (H) data. The 
inset (a) of Fig.l shows that the peak in J c at a given H 
manifests as an anomalous dip in R(t). As emphasized in 
Ref. 11, the location of each T p (H) is robust. The main 
panel of Fig. 1 shows portions of R(t) curves at different 
H and focusses attention onto the identification of t p (H) 
values. The inset panel (b) of Fig.l depicts the PE curve 
which passes through the t p (H) values for 120<H<380 
Oe. The shape of the PE curve clearly demonstrates its 
reentrant characteristics, with turnaround occurring at 
about 180 Oe, in excellent agreement with earlier results 
|l0f obtained via x'(T) measurements. The PE behavior 
in a given R(t) curve starts at a temperature, T p ; , where 
the resistance starts to decrease (see t p ; mark for H=360 
Oe curve). The data in the main panel of Fig. 1 imply 
that significant broadening of PE occurs as H values are 
decreased below 200 Oe, and eventually below 120 Oc, 
the PE becomes so shallow that t p (H) cannot be located 
precisely for H<100 Oe. 

A further confirmation of the reentrant characteristic 
of PE curve (in sample Y with T c (0)=7.17 K) is pro- 
vided by the results of dc magnetization hysteresis mea- 
surements, an example of which is shown in Fig. 2. In 
an isothermal scan at T=6.95 K, we expect to encounter 
both the lower and upper branches of the reentrant PE 
curve. The inset panel of Fig. 2 shows a portion of 
the M-H hysteresis data at T=6.95 K. The PE on the 
upper branch manifests as an anomalous opening of the 
hysteresis loop before reaching the H C2 value. According 
to Bean's Critical State Model |r7|| , magnetization hys- 
teresis AM(H) (difference between forward and reverse 
magnetization values at given H) is a measure of J C (H). 
To locate the existence of PE on the lower branch and 
to elucidate that the peak field H p on the lower branch 
of PE curve increases as T increases, we show in Fig. 
2 the variation of 4ttAM (cxJ c (H)) vs H at 6.9 K, 6.95 K 
and 7.0 K in the low field region. It is to be noted that for 
H< 60 Oe, AM(H) values at 6.95 K are lower than those 
at 6.9 K and this reflects the normal behavior in J C (H,T). 
However, a crossover in AM(H) curves occurs at about 
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70 Oe such that AM(H) values at 6.95 K exceed those 
at 6.9 K, thereby exemplifying the observation of anoma- 
lous behavior Q in J C (H,T) due to the presence of PE 
on the lower branch at 6.95 K. It may be further noted 
that AM(H) values at 7.0 K become larger than those at 
6.95 K above 130 Oe. Thus, at 7.0 K the (lower) peak ef- 
fect is seen to be encountered at a field value larger than 
that at 6.95 K, this establishes the central characteristic 
(i.e., dT p /dH>0) of the lower branch of the PE curve. 

Fig. 3 shows the temperature variation of the in-phase 
ac susceptibility (x' vs reduced temperature t) at a fre- 
quency of 211 Hz and in h ac of 1 Oe(r.m.s.) at some of 
the chosen dc fields Hdc (0 to 2 kOe) in three crystals 
X, Y and Z of 2H-NbSe2- These crystals were grown by 
the same vapor transport method |l9|] , but with start- 
ing materials of different purity such that the specimen 
X is the cleanest among the three and the specimen Z, 
the most disordered. The relative purity of the three 
crystals is reflected in the width of the superconducting 
transition in x'M response in zero field, which can be 
seen to progressively enhance in Figs. 3(a) to 3(c). In a 
x'(t) measurement at a given H, the PE is identified by 
an anomalous negative peak S|l3j. The peak tempera- 
tures T p (H) are independent of frequency and amplitude 
of ac field h ac in the parametric range of our experiment 
pof . Fig. 3(a) shows that in the most weakly pinned 
crystal X, the PE peaks are very sharp; their half widths 
are smaller than the width of superconducting transition 
in zero field. This fact supports the association of T p (H) 
with a possible phase transition, like, the phenomenon of 
change in spatial order of the vortex array across T p . 

Fig. 3(b) is based on the x'(t) response in crystal Y 
0- T P (H) values for H=300 Oe and H=100 Oe can be 
seen to be nearly the same, consistent with the reentrant 
nature of PE curve in crystal Y |l(| . However, the PE in 
H=200 Oe (data not shown) is significantly broader than 
that at H=300 Oe. Fig. 3(c) shows the x'(t) response in 
the most strongly pinned p3| crystal Z. In this sample, 
the PE peak is considerably broad even at H=1000 Oe 
as compared to the corresponding peaks in crystals X 
and Y (cf. Figs. 3(a) to 3(c)). The PE peak is barely 
discernible at H=300 Oe in crystal Z and below 200 Oe, 
the PE manifests only as a hump/kink in x'(t) curve. 

Fig. 4 summarizes the T p (H) data in three crystals X, 
Y and Z as H vs T p /T c (0) curves. For the sake of conve- 
nience and reference, we have also included in this figure 
the H Cl (t) and H c2 (t) curves obtained in crystal X [||. 
Thus, from Fig. 3 and Fig. 4, it is apparent that enhanced 
quenched random disorder (pinning) results in the fol- 
lowing effects : 

(i) For a given H (H>300 Oe), the PE peak occurs at a 
lower value of t (see, for instance, arrows marking T p (H) 
values in three crystals at H=300 Oe in Figs. 3(a) to 3(c)). 
This amounts to a "lowering" of the upper branch of PE 
curve with enhanced pinning (see Fig. 4). Also note that 
the PE broadens as the disorder increases. 

(ii) The reentrant lower branch of PE curve is clearly evi- 
dent only for crystal Y with intermediate level of disorder 



and not so for the other two crystals. The PE curve in 
crystal X shows a steep fall for 200 < H < 30 Oe and its 
turnaround characteristic presumably lies much below 30 
Oe; however, PE in x'(T) data is unobservable in crystal 
X for H<20 Oe. In sample Z, the PE becomes so broad 
at 300 Oe that the precise location of T p (H) below this 
field becomes somewhat ambiguous (see Fig. 3(c)). How- 
ever, it can be safely surmised that PE curve at lower 
fields (H< 500 Oe) moves away from the values that can 
be extrapolated from T p (H) vs H data at higher fields 
(H>500 Oe). Thus the effect of enhanced disorder (X 
to Z) is seen to result in "raising" to higher fields, the 
lower reentrant branch of the PE curve (in contrast to 
the upper branch). 

The results of Figs. 3 and 4 provide the following gen- 
eral conclusions: 

(i) For the upper branch (i.e., dT p /dH<0) of the PE 
curve, increasing pinning reduces the volume of Larkin 
domain of FLL, which then requires less thermal fluctu- 
ations to melt /amorphize it around t p . The quenched 
disorder and thermal fluctuations conspire together to 
destabilize the ordered phase and a "lowering" of the 
t p (H) curve is seen. 

(ii) For the lower branch (i.e., dT p /dH>0), increasing 
pinning reduces the volume of Larkin domain as before. 
Thus, one needs enhanced interaction (i.e., increase in 
H) to stabilize the ordered phase and the "raising" of the 
PE curve occurs. 

(iii) Increasing pinning has the general effect of masking 
the difference between a vortex solid and a vortex liquid, 
converting both into a "glassy" state. A clear distinction 
between the two phases is no longer possible, the trans- 
formation process is more gradual and thus a broadening 
of the PE results. 

(iv) As can be ascertained from the main panel of Fig.l, 
the onset (at t p z) of PE, which is apparently driven by 
quenched disorder, is also reentrant, just as the vortex 
melting transition is expected to be. Recent muon spin 
rotation studies |^l| on crystals of 2H-NbSe2 have pro- 
vided microscopic evidence in favour of a sudden change 
in spatial order of FLL at the onset of PE, i.e., at T p j. 

Finally, it is of interest to examine the scenario emerg- 
ing from Fig. 4 with the expectations of theoretical studies 
in the (H, T) region of our experiments. In the cleanest 
crystal X, where we observe only the upper branch of 
PE curve, which most closely marks the transformation 
of the ordered solid into the pinned "liquid" state and 
a narrow region separates this line from T C (H) line (cf. 
Figs. 6, 24 and 25 of Ref. 3 and our Fig. 4). In between 
T p (H)/H p (T) and T C (H)/H C2 (T) lines, the pinned liquid 
becomes unpinned at Hj rr (T) (see Fig. 2(a) for loca- 
tion of R lrr at 6.95 K). Fitting the higher field (H>200 
Oe) branch of PE curve to a standard quadratic Linde- 
mann relationship (see eqn. (1)) with Gi~3xl0 -4 and 
H C2 (0)«4.6 T Q, we get a reasonable value for Linde- 
mann number cl=0.15. A pronounced departure (in 
T p (H) at H<200 Oe) away from the upper portion ap- 
pears in accord with a recent theoretical work of Blatter 
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and Geshkenbein p2[ , who found a similar rapid drop 
in FLL melting curve (in the absence of pinning effects), 
away from the H C2 phase boundary, at low fields (i.e., 
above the turnaround feature of the melting boundary) . 
The reentrant lower branch of PE curve may thus be at 
still lower fields in crystal X, as proposed theoretically 
p2]| , and outside the detection capability of the present 
measurements. In crystal Y, with intermediate level of 
disorder, Ghosh et al [jlO| estimated Cl as about 0.17 from 
higher H portion of the T p (H) data. They also showed 
that if the Nelson-LcDoussal line Q| indeed marks the 
crossover from an interaction dominated to a disorder- 
dominated region, then the ratio of entanglement length 
Lb to the pinning length L c becomes approximately equal 
to 1 (L c /Lb~1) at B=30 Oe and t=0.975. L C /L E is given 
as 1, 



L c irK 2 ln(n) a$ 

~Le = ( V2 ( 2vrA(0) 



(J C )V2 (l- t )V3 



(G.jo) 1 / 2 



t 



(1) 



where various symbols have their usual meaning (sec cqn. 
6.47 of Ref.3). Below 30 Oe, the vortex array would be 
in disentangled liquid or glass state Eqn.l implies 

that the field at which crossover to glassy state can occur 
is cx Jc 1 / 2 - J c is larger in crystal Z by nearly a factor of 
50 as compared to that in crystal Y . This, coupled with 
the observation that in crystal Z the PE becomes diffi- 
cult to discern below t«0.96, yields a value of crossover 
field of ^400 Oe in sample Z, which demonstrates a very 
reasonable agreement with an observed value of 300 Oe 
in our experiment. 

In summary, we have demonstrated the variation of 
the PE at high temperatures and low fields, and espe- 
cially of its lower reentrant, branch is experimentally 
measurable. The results show not only the shifts in 
the relevant transitions/crossover lines separating the or- 
dered phase from the disordered ones, they also illustrate 
how the distinction between them is blurred as quenched 
disorder is varied. Some aspects of the results could 
be semi-quantitatively explained by the available theo- 
ries H| 22 2j| . Nevertheless, a more detailed theoretical 
analysis of the low field - high temperature regime, espe- 
cially on the liquid - dilute solid transformation process, 



is necessary. 
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FIG. 1. The inset panel (a) shows the temperature vari- 
ation of Resistance (R vs reduced temperature, t=T/T c (0)) 
in H dc =280 Oe (||c) in the crystal Y of 2H-NbSe 2 . The Peak 
Effect (PE) region can be identified by an anomalous dip and 
peak temperature t p is marked. The main panel shows por- 
tions of R vs t curves in the PE regions at some chosen dc 
fields in the range 140 Oe < H < 360 Oe. The marked t p (H) 
values elucidate that t p (H) increases as H decreases (dH/dt p 
< 0) from 360 Oe to 200 Oe and thereafter (for H< 200 Oe) 
dH/dt p > 0. The inset panel (b) focusses attention onto the 
reentrant nature of locus of t p (H). 
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FIG. 2. Isothermal magnetization hysterisis data [47rAM = 
4 7T (Mh-Mf)] for H dc ||c in the crystal Y of 2H-NbSe 2 at the 
temperatures indicated. The arrows identify the fields corre- 
sponding to onset of PE phenomenon on the lower (reentrant) 
branch of PE curve at temperatures of 6.95 K and 7.0 K. It 
is to be noted that the field at which (lower) PE occurs at 
7.0 K is larger than that at 6.95 K. The inset panel shows a 
portion of the M-H loop showing upper PE at 6.95 K. (Up- 
per) Peak field H p and irreversibility field H; rr at 6.95 K have 
been identified. 



FIG. 3. In phase AC (f=211 Hz, h ac = 1 Oe (r.m.s.)) sus- 
ceptibility (x') vs reduced temperature (T/T c (0)) in crystals 
X, Y and Z of 2H-NbSe2 in fixed dc fields (|| c) as indicated. 
The arrows mark the t p values at H dc =300 Oe in different 
crystals. 



FIG. 4. Magnetic phase diagram in three crystals of 
2H-NbSe2. PE curves and H C 2(T) curve correspond to 
T p (H)/T c (0) and T c (H)/T c (0) values obtained from x'OO 
data as in Fig. 3. For crystal Y, the two open squares ly- 
ing on the upper and lower branches of PE curve identify the 
peak field (H p ) and crossover field at 6.95 K as in Fig. 2. 
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Fig. 4 (S. S. Banerjee) 



